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Abstract 

We extend the validity of Dain's angular-momentum inequality to 
maximal, asymptotically flat, initial data sets on a simply connected 
manifold with several asymptotically flat ends which are invariant un- 
der a U(l) action and which admit a twist potential. 

1 Introduction 

In [3] Dain proved an upper bound for angular momentum in terms of the 
mass for a class of maximal, vacuum, initial data sets with a metric of the 
form 

g = e-2f^+2° i^dp^ + dz^) + /j^e"^^ [d^ + pBpdp + A.dzf , (1.1) 

where the functions are assumed to be (/^-independent. The existence of 
the global coordinate system (jl.ip has been justified for asymptotically flat 
axi-symmetric initial data sets on a simply connected manifold in the first 
paper of this series [2] . 

In this paper we extend the validity of Dain's inequality to all maximal, 
asymptotically flat, simply connected initial data sets (M, K) invariant 
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under a U(l) action, with several asymptotically flat ends and positive scalar 
curvature, and admitting a twist potential to as defined by (j2.6p below. 

In order to give a detailed statement, some preliminary remarks are 
in order. We choose an asymptotic region, say Mi and, following Dain, the 
remaining asymptotic regions are described by punctures on the z-axis in the 
(p, z) plane. The axial symmetry of the problem implies that the angular 
momentum Ji of each asymptotic region Mj is aligned along the rotation 
axis, and we shall write Ji for the relevant component of Jj. As is well 
known, iv is constant on each connected component of the punctured axis, 
and (cf., e.g., [9, Section 6]) Ji is proportional to the difference of the values 
of LO on the extreme segments of the axis, while for i ^ 1 the z-component 
Ji of the angular-momentum vector is proportional to the jump of uj at the 
i-th puncture. This implies 

N 

Ji = Y^Ji, (1.2) 

i=2 

so that Ji is determined by the remaining angular momenta. 

We denote by 87r/(J2, . . . , Jn) the numerical value of the action func- 
tional (j2.13p of the harmonic map, from = 0} to the two-dimensional 
hyperbolic space, constructed in Proposition 12.11 below. 

It follows from (jl.2p that for = 1 the angular momentum necessarily 
vanishes (so that no non-trivial inequality involving the angular momentum 
can be obtained in the = 1 case), while for A^ = 2 we have 



Ji = J2 and /(J2) = /(Ji) = 




We are ready now to present our version of Dain's inequality: 

Theorem 1.1 Let {M,g,K) be a three dimensional \J {1) -invariant initial 
data set, with positive matter density, on a simply connected manifold M 
which is the union of a compact set and of a finite number of asymptotic 
regions Mi, i = 1,...,N, N > 2, and with {g,K) — asymptotically fiat 
on each end in the sense of i2. 1\ }- JK^) with k > 6, together with 1^2.4^ . If 
(M, g, K) admits a twist potential satisfying \2.12\) . then the ADM mass 
mi of Ml satisfies 

mi > /(J2, . . . , Jn) ■ 
For N = 2 this is Dain's inequality mi > \J\Ji\. 
^This will necessarily be true in vacuum. 
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Conceivably the case of main interest is = 2 (already analyzed un- 
der different hypotheses by Dain [3]), since the function / is not known in 
general. (It would thus be of interest to obtain some lower estimates on / 
when > 3.) Now, it is far from clear how large is the class of metrics 
considered by Dain in [3], in view of several a priori restrictive conditions 
imposed there; the analysis in [2] and here shows that the inequality applies 
in substantial generality. 

A sketch of the argument seems appropriate: Following Dain, we show 
that the mass is bounded from below by the action of a map {U,uj), with 
values in the hyperbolic space, determined by the norm of the rotational 
Killing vector and by the twist potential. This map is singular on the rota- 
tion axis, with further distinct singularities on punctures on that axis, which 
correspond to the remaining asymptotically flat regions. One then wishes 
to show that the action of the map is bounded from below by the action of 
the extreme Kerr solution when = 2, or by the action of a harmonic map 
{U,oj) with singularities at the punctures resembling those of the extreme 
Kerr solution for N > 2. Such maps are constructed in Proposition 12. It the 
result is essentially due to Weinstein [10]. The key element of the remainder 
of the argument is a result of Hildebrandt, Kaul and Widman [6] that, on 
compact domains with smooth boundary, harmonic maps with negatively 
curved target space are minimizers of the action. Since we are working on a 
non-compact set, with maps satisfying singular boundary conditions, some 
work is needed to apply this result. We start by showing that the action is, 
roughly speaking, decreased by deforming {U,uj) to a map {Us,(l>s) with a 
singularity structure at the punctures resembling that of the extreme solu- 
tions. The maps {Urj,uJn) of Lemma 12. 51 are further deformations of (Us,ujs) 
which coincide with (U,ui) near the punctures, and for large r. Lemma 12.61 
introduces a final deformation which takes into account the fact that (U, cD) 
satisfies the harmonic map equations away from the rotation axis only. 

2 Angular-momentum inequalities 

We will consider Riemannian manifolds (M, g) that are asymptotically flat, 
in the usual sense that there exists a region Mcxt C M diffeomorphic to 
\ B{R), where B{R) is a coordinate ball of radius R, such that in local 
coordinates on M^xt obtained from \B{R) the metric satisfies the fall-off 
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conditions, for some k > 1, 

9ij - 6ij = Ofc(r-i/2) ^ (2.1) 

^kg^J e L2(Mext) , (2.2) 

R'jM = o(r-5/2) , (2.3) 

where we write / = Ok{r^) if / satisfies 

dk,...dkj = o{r^'') , 0<£<k. 

Let (M, g, K) be a general relativistic, not necessarily vacuum, initial 
data set, with 27r-periodic Killing vector r/. We impose asymptotic flatness 
on g as before and, in each asymptotic region, we assume the following 
asymptotic decay of K, for large r 

\K\g = 0{r-^) , A > 5/2 . (2.4) 

We will further suppose that A < 3, as faster decay would necessarily lead 
to zero angular momentum; A = 3 is the decay rate corresponding to the 
Kerr family of solutions. Note that ()2.4p enforces the vanishing of the ADM 
momentum of the initial data set. 

We also assume that the initial data set is maximal, iTgK = 0, and that 
the Einstein constraint equations hold with matter density ^ satisfying a 
positivity condition, 

^^^R=\K\] + lQ^T^x>\K\], (2.5) 

A key restrictive hypothesis in what follows is the existence of a twist 
potential to: 

eijkK^ ir]^r]Ux' = duj . (2.6) 

As discussed in [3] this holds e.g., for vacuum initial data sets on simply 
connected manifolds. 

It has been shown in [2] that, under the hypotheses of Theorem ll.H 
there exist coordinates in which the metric can globally written in the form 
(jl.ip (compare [1,5]). Consider an orthonormal frame such that 63 is 
proportional to 77, let 9^ be the dual co-frame; for definiteness we take 

Qi ^ e-^+'^dp , = e-^+^dz , 9^ = pe~^ {dip + pBpdp + A,dz) . 

By 

duj = eijkK^ trl'v^dx^ 

= e{dA,dB,d^)K{dx^ ,d^)dx^ 
= g{r^,v)e{ea,eb,e3)K{9\e3)9\ 
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for some (/^-independent function uj. Here, as before, the upper case indices 
A, B = 1,2 correspond to the coordinates {p, z), while the lower case indices 
a,b = 1,2 are frame indices. Thus, writing Kh^ for K{eb,e3) = K(0^,ez), 
we have 

p'^e-^^{K23e^ - KisO^) = dpuj dp + d.ujdz ; 



equivalently 



„3(7-a 



K 
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-d:,U) , K; 
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p3C/-Q 



SO that 



23 J 



(2.7) 
(2.8) 



In [2] (compare [1,3,5]) it has been shown that 



m 



1 

16^ 



(3)^ + ip2g-4a+2C/ (^^^^^ _ ^^^^^)2 J g2{a-t/)^3. 



+ 



(2.9) 



Inserting (|2.8p into (j2.9p we obtain 

771 > 



> 



1 
1 

8^ 



(3)i? + 2(Z?t/)2je2(-^)d3^ 



(i?[/)^ + — (Du;) 
P 



(2.10) 



It immediately follows from (j2.7[) that the twist potential uj is constant 
on each connected component j = 1, . . . , A^, of the axis £/ = {p = 0}; 
in this section we assume that N > 2. We set 



LOi := (J 



(2.11) 



As in [2, 3], in the coordinate system of p.ip each asymptotically flat end, 
except the chosen one, is represented by a point Si lying on the symmetry 
axis £/. 

To gain some insight into the problem at hand the following comments 
are in order: Roughly speaking, asymptotic flatness implies that the twist 
potential lo approaches some smooth function of the angles cjj, typically 
different in each asymptotic region, as one recedes to infinity there: 



iO — UJi = o(l) . 



(2.12) 



5 



The exact form of uJi is actually irrelevant for our purpose^; the essential 
point is that, after mapping infinity to a neighborhood of a puncture Si, dco 
behaves as l/r^ in the local coordinates there. This is at the origin of Dain's 
mass inequality: Indeed, the first term in ()2.10p is minimized by C/ = 0. 
However, a constant U would lead to an infinite mass integral because of 
the second term in (j2.10p . So the second term, plus the requirement that the 
integral converges, forces U to approach minus infinity as one approaches 
each puncture, enforcing thus a non-zero lower bound on the mass. 

Note that the explosion rate is faster for non-extreme solutions, compare 
(|2.14p and (j2.15p . giving a larger contribution from the first term, as com- 
pared to an extreme solution. However, there is a tension between the two 
terms, as the decrease of energy of the first term appears to be comparable 
to the increase of the second, which makes the comparison delicate. In fact, 
it is not at all clear a priori which term will win in the balance. In the 
argument below (closely related to, but not identical to Dain's [3]) we follow 
Dain's insight that cylindrical ends have less energy than asymptotically fiat 
ones: 

Proof of Theorem ll.lt If the mass is infinite there is nothing to prove, 
otherwise by ()2.10p we need to find a lower bound on 

I ■- j y^DUf + ^{DLof^d?x . (2.13) 

In a coordinate system where a Kelvin inversion has been performed in 
all remaining asymptotically fiat regions, as in [2], near each puncture the 
function U has, for small r^, the asymptotic behavior as in [2, Theorem 2.9], 

[/ = 21nri + 0(l) , dU = 2d{\nri) + 0{l) , (2.14) 

while ui approaches a non-trivial angle-dependent function uji as rj ap- 
proaches zero. 

We note the following result; when N = 2 the solution {U,uj) below is 
the pair of functions (U,uj) corresponding to an extreme Kerr metric with 
angular momentum along the z-axis equal to {uj2 — u>i)/8: 

Proposition 2.1 For any set of aligned punctures Si and of axis values ujj 
there exists a solution {U,oj) of the variational equations associated with the 
action ^2.13\) . with finite value of I, satisfying Ii2.11\) . with the asymptotic 
behavior near each puncture 

U = lnri + 0{l) , (2.15) 



^ An explicit form of uj for Kerr metrics can be found in Appendix 1X1 
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lD = 0(1) 



(2.16) 



Remark 2.2 It would be of interest to study in detail the regularity of the 
solution near the punctures, compare [8] for some related results. 

Remark 2.3 As discussed in Appendix[Cl an identical proof gives existence 
of harmonic maps with any prescribed number of non-degenerate (as in [9]) 
and degenerate horizons, with several asymptotically flat regions. We also 
prove there uniqueness of solutions. 

Proof: This is a rather straightforward consequence of the results in [10], 
we give some details to justify the estimates. Let the reference map (C/, u)) 
be any map from \ {oj} such that 

1. id takes the prescribed values uoj on the relevant connected components 
of the punctured axis of symmetry; 

2. {U,io) coincides, for all large values of + z^, with the extreme Kerr 
solution with the 2;-component of the angular-momentum vector equal 

to {ujn — '^i)/8; 

3. in a small neighborhood of the i-th puncture, after perhaps a constant 
shift in id, (U, ui) coincides with the extreme Kerr solution with angular 
momentum (wj+i — uJi)/8 near its cylindrical end; 

4. from the explicit form ()A.2p of the twist potential u> = ujkcit for the 
Kerr metrics one has, for small p, uniformly in z, away from the plane 
{z = 0} (see dHU-dAll) 

|w-Wi|<C^ and \Dlo\s<C^, (2.17) 

when (D = WKerr and n = r; one can therefore also arrange that (I2.17P 
be satisfied by the reference map ui away from the planes {z = aj}. 

5. To make things precise, near the axis p = and away from small 
neighborhoods of the punctures we let {jJ,Cd) be defined by the usual 
convex linear combination of two solutions using a smooth cut-off func- 
tion which depends only upon z for small p. 

For i large let {Ui,oji) be the map which coincides with ([/, u)) for /? < 1/i 
and for r > i, and such that {x := U — Inp, u)) solves the harmonic map 
equations, with target manifold metric 

b:=dx^ + e'^''duj\ (2.18) 
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away from the union of those last two sets; such a map exists by, e.g., [6]. 

By construction the tension map associated to {U — lnp,(D), as defined 
in [10], has compact support on M'^ \ {oj}, and is uniformly bounded in the 
norm defined by the metric h. Indeed, the last property is clear away from 
the axis. In the interpolation region near the axis the map (x, uo) is of the 
form 

x = -\np + aQ{z) + a2{z)p^ + 0{p'^) , uj = uJi + I3i{z)p'^ + 0{p^) , (2.19) 

for some smooth functions aoC-^), ol2{z)^ Pa{z)-, with the obvious associated 
behavior of the derivatives. The non-vanishing Christoffel symbols of the 
metric h are 

uju} 1 xu) U)X ^ ■ 

This leads to the following formula for the norm squared of the tension, 

\T\l = (Ax - 2e^^|L'w|2)2 + e^^(Aw + ADx ■ Duof , 

where A is the flat Laplacian on R^, with the scalar product, and norm of 
-D, taken with respect to the flat metric on M^. A uniform bound on |T|b 
readily follows from (j2.19p . 

Note that {Ui,u)i) has finite action / which is smaller than or equal to 
the action of (U,uj), as the action of {Ui,u)i) is strictly smaller than that of 
(C/,a)) on the region where they differ by [6]. 

As outlined in [10, Section 3], an appropriately chosen diagonal subse- 
quence of the sequence (C/j,a)j) converges uniformly on compact subsets of 

\ {oj} to the desired harmonic map {JJ — In p, uj), with {U — In p, to) lying 
a 6-finite distance from (U — ln/3,a)). The estimate (I2.15P follows. 

The action / of the limit is smaller than or equal to that of {U,lD) by 
Fatou's Lemma, in particular it is finite. □ 

The arguments in [10], together with elementary scaling in coordinate 
balls of radius p/2 centered at {p,z), show that there is a uniform gradient 
estimate 

\d{U- In p)\'^ + e^^^-^''P^\dQ\'^ = |(i([/-lnp)|2 + ^|dcD|2 < Cp-^ . (2.20) 

An identical estimate (with possibly a different constant, independent of i) 
holds for the approximating sequence, which implies that doj vanishes on 
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the punctured axis, and a) attains the desired values there. In fact, near 
from (j2.20p one obtains 

\dU\ < — , \doj\<^ . (2.21) 
p 

For further purposes we wih need a stronger estimate, which we prove in 
integral form. We consider small, non-overlapping balls near each puncture. 
Since all the functions are invariant under rotations around the rotation axis 
it suffices to work in a half-disc 

D^(C^b) := {x^>0, {x^f + (x2)2 < C^5} C 

of radius C\/5, with polar coordinates centered at (0, flj): {x^,x'^) = {psin6,ai+ 
pcosO). The reader is warned that the polar coordinate p here corresponds 
to in the applications that follow, and x^ here is p in the applications 
below; this explains the weight x^ in the measure in (I2.22p . 

Proposition 2.4 Let p > 1/2, and let uj he as in the proof of Proposi- 
tion \2.1[ and let 6 > be such that the half-disc D~^{^/5) centered at the 
origin contains only one puncture. There exists a constant Ci, indepen- 
dent of 5, such that for any positive measurable function g = ^dxl), where 
\x\ = \J {x^Y + [x^Y' , we have 

(lD-cj)" . 4^ x^dx^dx^ <Ci [ |L»(cD-lD)|2— g{\x\) x^dx^dx^' 



(2.22) 

Proof: Let, first, u be any function which vanishes near the axis x^ = 0, 
we claim that there exists a constant Ci, independent of 6, such that 

/ / xUx^dx^ <Ci [ \Du\^^^g{\x\)x^dx^dx^ . 

(2.23) 

In order to see that, we first prove that for a / —1, and for / G C^(0, tt), 
vanishing near zero and pi, 

r e-f{6 fde < r 9-+'f'{efde. (2.24) 

Jo [a + ly Jo 



9 



Indeed, 



a + 1 







2 



a + 1 



f{e)f{e)de 







The inequality (12.24p follows from the above. 

Next, for a < —2 (we will apply this with a = —2fi — 1, which then 
requires /i > —1/2), by using (j2.24p we obtain: 



< (^/2)" 



I e"^f{efde+ / [-K-eff^efdB 



/■TT /'TT 

f e''f{efde+ / e''f{^-efde 

JO 



< C'{a) r d^+\7r-eT+^f'{9f (useda + 2<0) 

JO 



(2.25) 
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Now, for X G D^{1) set /(x) = u{\^x), we have 
Jd+{V~S) sm^'^+^e \xy Jd+{i) sm^f'+^d \x\^ 

[p=o Je=o sm 



r f \ . \^^, g{p)p^ sin edpde 
Jp=o Je=o sm^^'+^ 9 

< Ci^6 C r \Dff-^g{p)p^ sin 9dpde 
Jp=o Je=o sm u 

= Ci [ \Du\^-^g{\x\)x^dx^dx\ 
Jd+(VS) sm^^'e 



< 



where in the step (*) we have used (|2.25p . 

Since uji coincides with u) for smah x^, we conclude that (|2.23|) holds 
with u replaced by — u). Passing to the limit z — > oo, (j2.22p follows. □ 

By an abuse of terminology, the couple {U,oj) constructed in Proposi- 
tion [211] will be referred to as an extreme Kerr solution; this is justified when 
there is only one singular puncture. 

Let {U,u}) be the functions U and uj given by Proposition 12.11 with the 
same value of cD on the axis as the map {U, to) under consideration, so that 

(w-ti)U = 0. (2.26) 

We will show that a lower bound on the action can be obtained by working 
in the class of f/'s of the form (j2.15p . For this let (5 > be small, we start 
by deforming {U,uj) to a pair (Lis, us) with the following properties: 

1. Away from balls centered at the punctures Oj of radius CVS, for an 
appropriate constant C, the new pair of functions {Us, (is) coincides 
with the original one {U,u>). 

2. For Tj < (5 the pair (Us,uJs) coincides with (Inr, d;), where tD is the 
function iv of Proposition 12.11 

3. The action / calculated for {tJs,(lis), which we denote by Is, is smaller 
than the action / calculated for the original solution, except perhaps 
for an error which tends to zero as 6 tends to zercH. 



^We expect Is to be strictly smaller than /, but this is apparent from our proof for 
small angular momenta only. 
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This can be done as follows: (j2.14p shows that for all < (5 < 1 small 
enough the equation 

U{ri,e) = \n5 
has a solution rj = fi{6, 5) ~ ^fS satisfying 

<5 < ri < e := C^fS , 

for a large constant C. We let Us to be equal to U away from a collection 
of non-overlapping balls centered at the punctures, where we set 

{Inrj, r-i < (5 ; 

In 5, 5<n<h{e,5)- (2.27) 
U{n,e), n>h{e,8). 

Then JJ^ is continuous, piecewise differentiable, hence in H^^^. Now, 

\DUs\^ = I \DUs\''+ I \DUs\^+ [ \DUs\^ 

JUi{0<ri<S}^ V ' JUi{5<ri<fi} ''r\i{h<ri} 





47r{N -1)6+ I \DU\'^ 



(recall that there are asymptotically flat ends, hence — 1 punctures). 
On the other hand 



f = f \DU\'^+ [ \DU\'^+ [ 



{0<ri<<5} Ju,{5<ri<f,} Jr\i{fi<ri} 

> I \DU\^ + / \DU\^ 



> 



U,{0<ri<<5} — V — ■' Jr\i{ri<ri} 



I rr2(4 + o(l))+ / \DU\ 

IQtt {N - 1)5 + o{5) + I \DU\'^ 



= 127r(Af-l)5 + o((5)+ / |Z)i/5p • (2.28) 

It clearly follows for all 5 small enough that the first term in / will be 
decreased when U is replaced by f/^. 
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It remains to check that the possible increase of the second term in / can 
be controhed uniformly in 5. For this we need to understand the behavior 
of Lv near the axis. It is convenient to rewrite (12.71) as 



d,uj = -e-3^+V2^i3 , dpu = e-^^+V^^^^^ (2.29) 



Condition (j2.4p implies that there exists a constant C such that in each 
asymptotically flat region we have 

\Du;\s < Cp^r^^ . (2.30) 

Let {x"^) = (yO, z) be the symmetry-adapted coordinates which extend to 
infinity in the i'th asymptotic region. Performing an inversion (compare [2]) 

/-af = x^/|x|2, (2.31) 

near each puncture a, = (0, Cj) = {af) we obtain 

\Duj\s < Cp^r^-^ . (2.32) 

This shows that 

'—\Du:\' = iDuf + '—\D^f 

P JvJi{0<ri<CV5} P Jn,{n>cV5} P 



in,{r,>Cv^} P 

0(6^)+ / Vl^^l' • (2-33) 

Jni{r,>CV5} P 



Next, let Cos be equal to uj away from a collection of non-overlapping 
balls centered at the punctures, while in those balls we set 



uj{ri,9), ri > C\fl ; 

u;(r,,^)^4Sr + '^(^-^)r7^' ^<n<C^b; 
uj{ri, 9), n < 6. 



We have 



P JUi{0<ri<5} P 



JUi{5<ri<cVS} P Jr\i{ri>cV5} P 
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The first term goes to zero as 6 goes to zero because {U , ui) has finite action. 
We claim that A goes to zero as 5 goes to zero as well, this requires some 
work. For 5 < < we rewrite ojs as 

Cbs{r„6) = (a; - Co){r„6)^^^^ + {Co - o)iri,0)Pj^ + 0{u,9) , 

where ui is as in the proof of Proposition 12. 1[ By (j2.27p and (12.14p 

e [n, f^j - I ^4C/(n,^^) < c'rf < C"6\ h{0, S) < n < CVS. ^ "^^ 
Hence, for some constant C2, 

C^'A < [ ^(\Du;\^ + \DCo\^ + \DCo\^ 

Ju.,{5<ri<C^/S} P ^ 



+ 



6^ {lo - Lof 



iu,{5<n<Cv^} P rf 
The integral involving Dui goes to zero as 6 goes to zero because 

-\Dcu\^<C3 —\m^, (2.36) 

Ui{5<ri<CV5} P JUi{5<ri<C^/5} P 

while {U,u)) has finite action; similarly for that involving Du). The integral 
involving Duj goes to zero by direct estimation using (j2.32p . Next, using 
Proposition 12.41 with g = 1 and /i = 2, we can write 



5^ (cD - cj)2 f 



\Ji{5<ri<CV5} P'^ rf Ju,{S<ri<CVS} P^ ' 

4 /■ 4 



The right-hand-side is integrable over the set Ui{rj < e} as in (j2.36p . and 
thus goes to zero as 6 does, hence also the left-hand-side. 

Consider, finally, the integral in the second line of (j2.35p . It is convenient 
to split the integration region into two, according to whether or not |z — aj| < 
p. In the region 

n -.= {5 <ri<C^5 , \z-a,\ < p] , 
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the function p is equivalent to rj, while both oj and a; are bounded there. 
This gives the straightforward estimate 



In the region 



— 2 = ^ • 



^1 ■■= {6 <ri< CVS , \z - Oil > p} 



the function — aj| is equivalent to r^. Both uj and u) satisfy (j2.32p . By 
integration along rays within the planes z = const from each connected 
component of the axis we obtain in Yi, for ri < e small enough and 
z 7^ Oj, 

\oj - Cj\ < C p-^\z - ai\^^'^ . (2.37) 

Integration over Yi gives 

which goes to zero by our hypothesis that A > 5/2. 
Summarizing all this, 

I>Is- o(l) , (2.38) 

where o(l) goes to zero as 5 does. Smoothing out (lis, Cos) at the corners, 
we can further assume that {Us, Cos) is smooth without affecting ()2.38|) . 

To continue, we show that Is > /, where / is the action corresponding 
to the map of Proposition 12. li In order to prove this, let e > be such that 
the balls of radius e centered at the punctures do not overlap, and note that 
there exists a large constant C such that both {Us, Cos) and {U,Co) belong to 
the class .^^a of maps {U, uj) defined as follows: 

1. / [\DU\^ + ^\Dio\^] < C; 

2. For < Tj < e we have \U — \nri\ < C; 

3. Lo coincides with Co near the punctures; 

A. LO = LOj on 

It also follows from [2] that the map {U,lo) of Theorem 11.11 satisfies the 
following 

— 1/2 

1. For < Tj < e we have \DU\s = o{r- ); 
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2. for < /O < e, we have, 

i 

for some 5/2 < A < 3 and for each i, where the first term r}~^ accounts 

— A 

for small behavior, while the second one accounts for large r's; 

3. We have U = o(r"V2) j^^f/j _ Q^j, 3/2^ ^ tends to infinity. 

(this is also true for ([/, Co) when N = 2, and is expected to be true without 
this restriction on A^, but such an estimate has not been established so far). 

Dropping the subscript {Us,ujs) ^or notational simplicity, we wish to show 
that the corresponding action / is larger than or equal to that of (U,u)). In 
order to see that, let > and let cpr^ G C°°(M'^) be any family of functions 
satisfying 

• < ipr^ < I; 

• ifrj = for < rj < r//2 and for > 2/ry; 

• iprj = 1 on the set Wr^, where 

Wx := {r < 1/A} {r, > A} ; 

• I-D'/'tjI < C/rj for i]/2 < ri < rj; and 

• l-Dv^ryl < Ct] for 1/r] <r< 2/r]. 

Set 

Un = ipr^tj + (1 - ^r))U , 
(jj^ = LPrjCo + (1 - ^r^Cb . 

Note that, by definition of cD, we have ojrj = Co near the punctures for r] small 
enough. Using again / denote the value of the action / for (f7,cj), we claim 
that the action of (C/^,a;^) satisfies 

Lemma 2.5 lim^_»o Ir^ = I- 
Proof: Indeed, we have 
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= / \Dtj\'^+ / \DUr^\'^ + 

JU,{0<r,<V2} •^Ui{r?/2<ri<r;} 



/ // 

+ / \DU\'^+ [ \DUX+ I \DU\^ 

r\i{r,<Ti , r<l/»?} J{^h<r<'2h} J{2/v<r} 



III IV V 



The integrals / and V converge to zero by the dominated convergence the- 
orem, while /// converges to the integral over of |L>f7p by, e.g., the 
monotone convergence theorem. 

The term IV can be handled as follows: 



[ \DU^\'^ = [ \{U - U)D^r^ + ip^Dfj + {1 - ipr^)DU\'^ 

J {l/v<r<2/v} J{^/ri<r<2/ri} 

< 3/ (\{U - U)D^^\^ + \DU\^ + IDUI"^ 

J {llri<r<2ln\ ^ 



' {l/ri<r<2/r)} J {l/v<r<2/r]} 

f 

'{l/r}<r<2/-n} 

< 3 /" (C{U -Ufr^^ + \DU\^ + \DU\'^ 

J{l/ri<r<2/ri} ^ 

The second and third term go to zero by the Lebesgue dominated conver- 
gence theorem. Letting {U,lD) be as in the proof of Proposition 12.11 for 
small T] the first term can be estimated as follows: 

r-\U-tjf = [ r-\U-Uf 

{l/r]<r<2/ri} J {I / ri<r<2 / r)} 

< 2 [ r^^ilU^ + 2U^ + {tj -U)'' 

J{l/v<r<2/ri} ^ 

The first two integrals tend to zero by direct estimations. As U — U is 
the limit of compactly supported functions the weighted Poincare inequality 
applies to the third term, implying that the function r~^(C7 — f/)^ is in L^. 
The vanishing of the limit of IV as r] goes to zero follows now from the 
Lebesgue dominated convergence theorem. 
The analysis of // is identical. 

A similar analysis applies to the remaining integral in The only 
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delicate term is 



4 \-^^v\ - / \-^^ri 
{l/ri<r<2/'ri} P J {l/v<'r<'i/ri} 



f c 

< / —\{UJ - Oj)Diprj + iprjDuJ + {1 - ipr^)DuJ 
J {l/ri<r<2/ri} P 

< 3/ ^ {\{uj - u;)D^r,\'^ + \Duj\'^ + IDCjI"^) 

J{l/n<r<2/r]} P 

^ ~.\2-2 , I .|2 , I T-,r.\2\ 



< 3/ — (C(w - + |L>wp + |L>tD 

J{l/ri<r<2/'ri} P 

The last two terms go to zero as before. The first can be estimated as 

{l/r]<r<2/ri} P ^ •^{l/'?<r<2/r)} P r \ / 

The first term goes to zero by direct estimations. The weighted Poincare 
inequality (j2.22p with /x = 1 and g{r) = applies to the last term, giving 



{i/v<r<2/ii} P^r'^ J{^/v<r<2/-n} ^''^ sin^ - 

< / \\D{u-Co)\' 

J{l/'ri<r<2/ri} P 

r /p4f7 „4C7 

< C r—\DCo\'+'—\DCo\' 

J{l/rt<r<2/r)} ^ P P 

and the right-hand-side goes to zero in the limit. 

□ 

The next step is to prove 
Lemma 2.6 Irj > I for all rj small enough. 

Proof: For all rj small the maps {Urj,uJrf) and {U,oj) coincide on balls of 
radius r//2 around the punctures, as well as on the complement of a ball 
of radius 2/7]. One would like to use a result of [6], that the action / 
is minimized by the solution of the Dirichlet problem, which is expected 
to be {U,uj); however, that result does not apply directly because of the 
singularity of the equations at the axis p = 0; moreover, we are working in 
an unbounded domain. To take care of that, for e < 1 let 



0, p < e; 
ln£ 

In^ 

1, ' p>^. 
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Set 

C/^,£ = (pJJj^ + (1 - '^e)U , UJrj,e = ^e^^ri + (1 " ^e)^ ■ 

Let Ir^^e denote the action of {Urj^e-,^r},t)- We claim that 



f r e^c/^.E -1 

/ {DU^^.f + {Duj^^.f d^x^.^oO. (2.39) 

J{p<^e} L /5 J 



Equivalently, 

Ir,,e ^e-.0 Ir, ■ (2.40) 

In order to see this, note that, for e < 77/2, the integrand of (j2.39p is non- 
zero only away from balls of radius ri/2 centered at the punctures, with 
moreover r < 2/r]. Next, the integral over the set {0 < p < e} approaches 
zero as e tends to zero by the Lebesgue dominated convergence theorem. So 
it remains to consider the integral over 

■.= {e<p<V~e}n{r< 2/77} n, {u > r?/2} , 

which can be handled as follows: 

{DU^,ef = [ ({Ur,-U)Dip,+ip,DUr, + {l-iPe)Du]^ 



<C{p\ lne|)-i 



< 3/ ( +\DUr,\'^ + \DU\^ 



p2|Jj^g|2 ^ 



1 



In e 



0<p<e ,r<2/ri 



< 0{-^)+3l {\DUr,\^ + \DW^ 



The integral in the last line goes to zero by the dominated convergence 
theorem. 

The analysis of the term containing the derivatives of u;,,^e is similar, 
using (|2.32|) . (|2.37|) . and Proposition 12.41 with p = 4 and g{r) = r^: 

f g4C/j7,e f Ciri] / \ 2 

/ -^{Dujr,,er < / -^[{iOr^-Cj)D0, + ^,D{LOr^-C[;) + DC[;j 

^ sm t^l lne| V / sm V 

-^r,,. sm 



< o{e) + CC'{r]) [ ( \Duj\^ + \DLd\'^ + \DiJj\'^) ^,^oO . 
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This ends the proof of (j2.39p . 

Now, w^^e) coincides with {U,uj) on the set {p < e}, and on balls of 
radius i]/2 around the punctures, and on the complement of a ball of radius 
2/r]. Further, after shifting U by Inp, the variational equations associated 
with the action / are the harmonic map equations, with target space — 
the two-dimensional hyperbolic space. Hence the target manifold satisfies 
the convexity conditions of [6] (see Remark (i), p. 5 there). We can thus 
conclude from [6] that action minimizers with Dirichlet boundary conditions 
exist, are smooth, and satisfy the variational equations. It is also well known 
that solutions of the Dirichlet boundary value problem are unique when 
the target manifold has negative sectional curvatures. All this implies that 
(C/,a;), with its own boundary data, minimizes the action integral over the 
set 

{p>e}n{r< 2/r]} n^ {r^ > r]/2} . 



Hence 



d^x > 



DU) {Du 



~\2 



d^x . 



By the monotone convergence theorem we have 



/ 

■J{p>e} 



DU] +^—{Duf 
P 



d-'x 



SO that 



= lim 
> lim 

= i . 



e 



DU) +-^{D^f 



d^x 



d-^x 



□ 



Returning to the proof of Theorem 11.11 Lemmata 12.51 and 12.61 applied 
to (f7,cj) = {JJsTibs) give Is > /. Passing to the limit (5 — > we obtain, 
by (|2.38p . I > I. In the case of two asymptotic regions one concludes by 
noting, following [3], that 



□ 
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A Kerr solutions 



The Kerr black holes provide an exphcit family of solutions of the singular 
harmonic map equations, as follows: In Boyer-Lindquist coordinates, which 
are denoted by (t, f, 0, 99), the metrics take the form 

A — sin^ 6 n Amaf sin^ 9 
g = at H dtdip + 

+ sm2 edip^ + —dr^ + T^dO^ . (A.l) 

Here 

E = f2 + cos^ 9 , IS. = r"^ + a? — 2mr = (f — r+)(f — r_) , 
and r+ < f < 00, where 

r-t = m lb (m — a ) 2 . 

The function oj reads [4] 

^, o „ ma^ cos 0sin^0 , 
w = J(cos^6' - 3cos6') . (A.2) 



The constant m is of course the total mass and a = J /m, where J is the 
angular momentum. Note that the leading order term in u) is uniquely 
determined by J. 

We relate now f and 9 to p and z. If \a\ < m let r_|_ = m + ^frr^^^^a?' be 
the largest root of A, and let r+ = otherwise. For 

f > r+ , 

so that A > 0, define a new radial coordinate r by 

r = ^ (^f - m + \/A^ ; (A.3) 

this has been tailored [4] so that after setting 

p = rsin0, z = rcos0, (^-4) 
the space-part of the space-time metric takes the form 

g = e-2C>+2« (^p2 ^ ^^2^ ^ p2g-2t/ ^ ^^^^^ ^ A^dz)^ . (A.5) 
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We have 



— 



f = r + m -\ . (A. 6) 

4r 

We emphasize that while those coordinates bring the metric to the form 
(jA.Sp , famiUar in the context of the reduction of the stationary axi-symmetric 
vacuum Einstein equations to a harmonic map problem, the coordinate p 
in ()A.4p is not the area coordinate needed for that reductiorS except when 
m = a. 

To analyze the behavior near r = we have to distinguish between the 
extreme and non-extreme cases. Let us first assume that m? 7^ a?, then we 
have 



U 





1-^ 


(-) - '° 






2 



+ 0(r). (A.7) 



On the other hand, in the extreme case m = a it holds that 

U = In (—) + - In (1 + cos^ 6) + 0(r). (A.8) 
\2m/ 2 ^ ^ 

Furthermore, for m = iba, for small r, 

a; = TW +0(r). (A.9) 

1 + cos^ 

We will also need derivative estimates for u: from ()A.2p we obtain the 
uniform estimates 

|5,a;|<Csin3e = c4, \dM < C^-^ = C f , (A.IO) 
so that, for r < 1, 



,3 



|I)u;|5 = + r-2(a,c^)2 < . (A.ll) 



r4 

We claim that away from the plane z = we have the uniform estimate 

\oj-uJi\ < . (A.12) 

where loi = —2J is the value appropriate in the half-space {z > 0}, while 
LV2 = 2J is the one which should be used for ^ < 0. This is clear for the last 
term in ()A.2P : for the first one this follows immediately from the following 
identity, obtained by expanding cos^ 9 = {cos 6 — cos(n7r) -|- cos(n7r))^ with 
n = 0,1: 

cos^ 9 — 3 cos 9 = —2 cos(n7r) -|- (cos 9 — cos(?7-7r))2(2 cos(n7r) -|- cos 9) , 

''The correct {p, z) coordinates for the harmonic map reduction are p = v^sin6', 
z = {f — m) cos 0. 
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B Asymptotic equations near a degenerate hori- 
zon 



The maps {U,(jj) given by Proposition 12 . 1 1 are axially-symmetric (i.e., invari- 
ant under rotations around the z-axis) solutions of the variational equations 
for the action 



1 

8^ 



2 p4f/ 

DU] +—-{Duf d^x. (B.l) 

p4 



where D is the usual gradient of the flat Euclidean metric on M'^. Thus the 
equations read 

A^7 = 2^(L>w)2 , (B.2) 




= . (B.3) 

\ P' ) 
One expects that 

[/ = In r + ?7(0) + 0(r) , cD = ^(0) + 0(r) , 

where the remainder terms are bounded under r- or ^-differentiation^ In- 
serting this into the equations ()B.2|) - ()B.3|) it must hold 

^a,(sin0a,f7) = -l + 2sin-4 0e4f>(a,cl;)2 , (B.4) 
de (sin-3 Oe'^^de^^ = . (B.5) 

Elementary ODE theory shows that these equations admit a four parame- 
ter family of solutions, some of which might fail to be regular everywhere. 
Note that at this stage the question of differentiality of the solutions at the 
axis is open, and any such regularity in our context would need a careful 
justification. (Recall, however, that for the solutions we are considering the 
functions U and uj are bounded.) 

Solving (jB.SP for uj one finds that there exists a constant c such that 

-^deisinO dell) = -1 + 2c^ sin^ 9e-^^ . (B.6) 



^In fact, assuming the weaker condition that any error term, say e(r, 0), has the property 
that de€ = 0{r), dr€ = 0(1), d^e = 0(l/r), 9|e = 0(r), leads to error terms 0(r) in 



(|R4| - (|R5| . and to errors 0(l/r) in (fR2ll - ((R3| when ((R4)) - ((R5| hold 
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All solutions with c = have a constant uj and, for some constants a, /? G M, 
U = Q + (1 - /3)lnsin6l + /31n(l - cos6l) , 

so that 



U = Inp + a + /?( ln(l — cos 



In sin 9 



This is not of the form we are looking for as the correction U to In r is never 
bounded. 

Next, the addition of a constant to U can be used to rescale the constant 
c to an arbitrary value; e.g., one can choose 2c^ = 1. 

Finally, let a,/3 S M; one readily checks that the couples {U,u>) given by 



U = a + ^ln(l + cos2 6') 



u> = ibe" 
satisfy (lR4l) - (lR5ll . Note that 

douj = 



-2a 



COS 9 



1 + cos2 9 



+ /3 



(B.7) 
(B.8) 



e ^" sin^ 9 
(1 + cos2 0)2 



which exhibits clearly the high order of vanishing of di°o at the axis. This 
agrees with the functions arising from ()A.8p - (jA.9p if 



a 



ln(2m) 



According to Jezierski [7], the general solution of ()B.6P with c 7^ can 
be parameterized by two real constants (3 and 7, with (3^ ^ 0, as follows: 



[/ = - In 
2 



(1 - cos 9f ± 7^(1 + cos 9f] sm^-P 9 



Pi 



(B.9) 



In fact, it is straightforward though tedious to show that these functions do 
indeed solve ()B.6p . Further, within this family we have 



1 



In 



c(l±7^ 



deU{ 



2(1 ±7^) 



One checks that this can always be solved for (3 and 7 in terms of 

and d0U{^), showing that all solutions of (|B.6p which are bounded near tt/2 
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are given by ()B.9|) . Note that the case /? < in ()B.9P can be reduced to 
/3 > by simple redefinitions. So, assuming /? > 0, one checks that the only 
solutions which are uniformly bounded are the ones with /? = 2, 7 > 0, and 
with the plus sign chosen 



U=-ln 
2 



c((l-cose)^ + 7^(l + cos(6'))^ 
2^ 



(B.IO) 



One can then integrate ()B.5|) to obtain uj. We thus conclude that there 
exists a two-parameter family of bounded solutions of (|B.4p - ()B.5p . and that 
the extreme Kerr solutions provide only a one-parameter family thereof. 

If one approximates {U,Cd) by (Inr -|- U,Lb), then ()B.2p - ()B.3P will be 
satisfied up to terms 0(l/r). It should be clear to the reader how to push 
the expansion one order higher to obtain a bounded tension map, but the 
details of this calculation have no interest. 

In the notation of [9], the target space metric takes the form 



In this parameterization, (|A.8P - ()A.9p can be rewritten as 

sin^ 6 =f2 cos ( 



X,Y 



Am' 



1 + cos2 1 + cos2 e 



+ 0{r) . 



(B.ll) 



(B.12) 



In this context one could also look for solutions which depend only upon 
9] the equations then read 

dei^iTiOX-^deY) = , -^de (sin^^eX) = X-^UdeXf - {deVf] . 

sm d V / 

A solution is given by 



Thus Idyl 



X,yj = (^sin(0),cos(0) 
\dQY/r\ equals p/r"^, as in the scaling estimate 
\dY\ < C4 , 



(B.13) 



and not better. However, the solution ()B.13P does not behave like the 
extreme solutions we are looking for: U here equals Inr plus an angle- 
dependent correction as desired, but the latter blows-up badly at the axis. 
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C On uniqueness of harmonic maps associated to 
black holes 



It is expected that to a stationary "multi-black-hole" vacuum space-time 
one can associate a harmonic map which lies to a finite distance, in the 
hyperbolic target space, from a map with the following properties, modelled 
on a Kerr solution: 

1. There exists A^^^h ^ degenerate event horizons, which are represented 
by punctures {p = 0, z = bi), each of them labeled by a mass param- 
eter rrii > and angular momentum parameter Oj = =bmj, with the 
following behavior for small rj := -y^ p"^ + {z — bi)^, 

The twist potential u; is a bounded, angle-dependent function which 
jumps by 4Jj when crossing 6j from z < 6j to z > 6j, where Jj is the 
"angular momentum of the puncture" . 

2. There exists A'ndh ^ non-degenerate horizons, which are represented 
by bounded open intervals Jj C £/, with none of the previous bj's 
belonging to the union of the closures of the /j. The functions U — 
2lnp and uj extend smoothly across each interval /j , with the following 
behavior near the ends points, for some constant C: 

\U-^ ln(V/o2 + iz- CiY + z - Ci)| < C near(0,Q). (C.2) 

The function u; is assumed to be constant near the q'sIeI 

3. The functions U and uj are smooth across -s^ \ (Uj{6j} Uj /j), with lo 
locally constant there. 

As pointed out by Dain, and used in our work above, an alternative 
way of representing a non-degenerate Kerr black hole is provided by a map 
into the hyperbolic space, which is not harmonic, with a puncture on the 
symmetry axis corresponding to the second asymptotically flat region. This 
generalises naturally as follows: 

®For the Kerr solution the twist potential lo is of course not constant near the end 
points, but this simple condition is good enough for the purposes of Theorem IC.ll 
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4. There exists A'af + 1 asymptotically flat regions, for some A^af ^ 
0. A set of explicitly asymptotically Euclidean coordinates for the 
first asymptotic region is provided by x = {pcosip, psiiup, z), with 
taking large values. The remaining asymptotic regions are represented 
by punctures 6j = (0, 0, 6j) G .s/. If we set 



then we have the following asymptotic behavior near each of the punc- 
tures, which corresponds to a Kelvin inversion of asymptotically flat 
coordinates of a Kerr solution with mass nii and angular momentum 
parameter Oj, 



(see [2, Theorem 2.9]; compare Appendix |A]) . The twist potential lo 
jumps by 4Jj when crossing bi from z < bi to z > bi. 

The structure described in points 1-4 above will be referred to as the axis 
data. Thus, some of the punctures bi correspond to degenerate horizons, 
while the remaining ones correspond to further asymptotically flat regions. 
Defining the distance between two maps and ^2 as 



where the distance dj, is taken with respect to the hyperbolic metric (j2.18p . 
we have the following generalisation of Proposition I2.lt 

Theorem C.l For any set of axis data there exists a unique harmonic map 
$ : \ — > Jif"^ which lies a finite distance from a solution with the 
singularity structure above, such that uj = on £/ for large positive z. 

Remark C.2 There does not seem to be any obvious relationship between 
the harmonic maps here with A^af 7^ 1 and stationary vacuum black holes: 
we emphasise that the map $ corresponding to non- degenerate Kerr black 
holes is not harmonic in conformal coordinates {p, z) in which the second 
asymptotically flat region is represented by a puncture. 

Remark C.3 For A^af = A'dh = existence, and uniqueness under a sup- 
plementary H"^ condition, have been previously proved by Weinstein [10]. 
Similarly, uniqueness under again an additional condition for A^af = 
A^ndh = 0, A^dh = 1 has been proved by Dain [3]. 



^p2 + _ J,^)2 ^ 




(C.3) 



sup 4(^>i(p),$2(p)) 
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Proof: Existence can be established by repeating the proof of Proposi- 
tion 12.11 For uniqueness, a simple proof can be given as follows: Be- 
cause of the negative sectional curvature of the target, the distance function 
f{p) = <l>2(p)) > is subharmonic on \ .s/. The vanishing of / 

follows then from Proposition IC.4I below. □ 

Recall that denotes the 2:— axis. We have: 
Proposition C.4 Let f e C^{R^ \ £/) satisfy 

Af > in \ £/, in the distribution sense, (C4) 

0</<l, onR^\£/, (C.5) 



and 



hm f{x,y,z)=0. (C.6) 

{x,y,z)eR^\s^/,\{x,y,z)\—^oo 



Then 

f = 0, onR^X £/. 

Proof of Proposition 1. Given any e > 0, there exists, because of (\C.Q\i . 
some positive constant R > 0, such that 

f{x,y,z)<e, y {x,y,z) gR:'\.^, \{x,y,z)\>R. (C.7) 

For < 5 < i?, let 

Ds := {{x,y,z) \ \{x,y,z)\ < R,\{x,y)\ > S}. 

Define, on Dg, 

, logi\{x,y)\/R) 
'■='^ log{6/R) ■ 

Clearly 

g5> e on Ds- 
In particular, in view of (|C.7p . 

g5{x,y,z) > e > f{x,y,z), V {x,y,z) G dDsr\{{x,y, z) \ \{x,y,z)\ = R}. 
Using (jC.Sp . we also have 

95{x,y,z) = 1+e > f{x,y,z), V {x,y,z) G dDsn{{x,y,z) \ \{x,y)\ = 6}. 
Thus we have proved 

gs > f, on dDs. 
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Since gs is harmonic in Dg, and / is subharmonic in Ds, we have, in view 
of the above, 

95 > f, on Ds. 

Namely, for the e and R, 

' + ^^'^IK^^ > fix, y, z), y\{x,y,z)\<R,\ix,y)\>6,R>6> 0. 

Sending 5 to in the above leads to, for the e and R, 

e>f{x,y,z), y\{x,y,z)\<R,\{x,y)\>0. (C.8) 
This, together with ()C.7p . implies 

e > /, on \ i/. 

Sending e to zero leads to 

> /, on \ £/. 

We have thus proved 

/ = on \ 

and the proposition is established. □ 
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